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On  the  Applicability  of  Directed  Fluid  Jet c tc 
Newtonian  and  Non-Newtonian  Flows 


by 


P.  M.  Naghdi 


Abstract.  This  paper  elaborates  on  the  applicability  of  a direct  formulation 
of  the  theory  of  jets  to  Newtonian  and  non-Newtonian  fluid  flow  problems. 
Following  a general  discussion  of  the  nature  of  direct  theories,  we  record  the 
basic  equations  of  the  theory  of  directed  curves  for  any  finite  number  of 
directors.  Reference  is  then  made  to  some  recent  results  for  incompressible 
Newtonian  viscous  flows,  including  the  problem  of  jet  breakup.  The  major 
portion  of  the  paper  is  concerned  with  application  of  the  direct  approach  to 
an  incompressible  non- Newtonian  Poiseuille  flow  in  a circular  pipe.  The 
results  are  compared  to  those  of  the  three-dimensional  theory  and  are  found 
to  include  the  effect  of  "normal  force"  corresponding  to  the  "normal  stress" 
effect. 
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1.  Introduction.  General  background. 

The  objective  of  this  paper  is  to  indicate  the  manner  in  which  a mechanical 
theory  of  fluid  jets  constructed  by  a direct  approach  (rather  than  from  the 
three-dimensional  equations)  may  be  used  to  treat  certain  one-dimensional 
Newtonian  and  non-Newtonian  flow  problems  which  may  not  be  tractable  with  the 
use  of  exact  three-dimensional  equations.  This  direct  approach  for  fluid  jets 
is  based  on  one-dimensional  models  called  Cosserat  (or  directed)  curves  (defined 
in  Sec.  2),  which  are  curves  in  a Euclidean  3-space  endowed  with  additional 
structure  in  the  form  of  kinematical  vector  fields  called  directors . Clearly, 
if  full  three-dimensional  information  is  desired  regarding  the  motion,  deforma- 
tion and  distribution  of  stresses  of  the  continuum  under  study,  then  there  is 
no  point  in  developing  a one-dimensional  theory.  In  fact  the  aim  of  a one- 
dimensional theory  of  the  type  considered  here  is  to  provide  only  partial 
information  in  some  specific  sense:  for  example,  in  the  case  of  a fluid  jet  or 
fluid  flow  in  a pipe,  information  concerning  quantities  which  can  be  regarded 
as  representing  the  medium  response  effectively  confined  to  a neighborhood  of  a 
curve  as  a consequence  of  the  (three-dimensional)  motion  of  the  body,  or  the 
determination  of  certain  weighted  averages  of  quantities  resulting  from  the 
(three-dimensional)  motion  of  the  body. 

The  development  of  the  basic  theory  of  Cosserat  curves  is  exact  in  the 
sense  that  it  rests  on  (one-dimensional)  postulates  valid  for  nonlinear  behavior 
of  materials.  By  tne  nature  of  its  construction,  the  theory  necessarily  satisfies 
the  requirements  of  invariance  under  superposed  rigid  body  motions  that  arise 
from  physical  considerations  and,  of  course,  is  also  consistent  and  fully 
invariant  in  the  mathematical  sense.  Moreover,  the  development  by  the  direct 
approach  is  conceptually  simple  and  is  free  from  the  difficulties  involved  in 
the  approximations  usually  made  in  the  derivation  of  jet  theories  from  three- 
dimensional  equations. 


1. 


It  should  be  remarked  here  that  the  use  of  a direct  approach  based  on 

Cosserat  curves  in  formulating  one-dimensional  theories  does  not  mean  that  one 

ignores  the  nature  of  the  field  equations  in  the  three-dimensional  theory.  In 

fact,  some  of  the  developments  of  the  field  equations  by  direct  procedure  are 

materially  aided  or  are  influenced  by  available  information  which  can  be 

obtained  from  the  three-dimensional  theory.  For  example,  the  integrated  three- 

dimensional  equations  of  motion  establish  guidelines  for  a statement  of  the 

conservation  laws  of  the  one-dinensional  direct  theories,  and  also  provide  some 

insight  into  the  nature  of  inertia  terms  and  the  kinetic  energy  that  appear  in 
* 

the  latter  theories. 

Such  difficulties  as  are  associated  with  the  derivation  of  one-dimensional 
theories  from  the  three-dimensional  equations  most  often  arise  in  the  construction 
of  constitutive  equations,  and  it  is  in  fact  here  that  the  direct  approach 
offers  much  appeal.  While  an  approximation  to  the  constitutive  equations  in 
the  three-dimensional  theory  retains  the  constitutive  coefficients  which  have 
been  predetermined  within  the  scope  of  the  three-dimensional  theory,  the  use 
of  such  results  in  a one-dimensional  theory  may  in  general  lead  to  incorrect 
results.  In  this  connection,  it  should  be  observed  that  the  constitutive 
coefficients  of  the  direct  theory,  in  general,  may  involve  contributions  from 
both  the  material  properti- s of  the  three-dimensional  medium  and  the  looal 
geometry  of  the  bod.y  (here,  the  jet-like  body)?  The  procedure  employed  in  the 
direct  theory,  on  the  other  hand,  leaves  the  constitutive  coefficients  unspecified 
and  while  the  determination  of  these  coefficients  may  require  substantial  effort, 
they  can  eventually  either  be  related  to  those  of  the  exact  three-dimensional 
theory  or  else  be  determined  by  suitable  experiments. 

*See  also  the  remarks  following  equations  (2.15). 

^As  an  illustration,  see  for  example  Eqs.  (62)  of  [4]  in  which  the  constitutive 
equations  for  an  incompressible,  viscous,  elliptical  jet  depend  not  only  on  the 
shear  viscosity  buc  also  on  the  local  time -dependent  geometry  of  the  cross- 
section  of  the  jet. 
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i. 


As  defined  in  Appendix  A,  a jet  is  a three-dimensional  body  whose  boundary 
surface  has  special  features  and  to  this  extent  it  is  similar  to  a rod  although 
the  nature  of  the  specified  surface  (or  boondary)  conditions  in  the  two  bodies 
may  be  different.  Moreover,  the  kinematics  of  jets  and  rods  are  identical  and 
it  is  only  through  the  constitutive  equations  that  a distinction  appears  between 
rods  and  jets.  In  many  ways  the  development  of  the  theory  of  Cosserat  curves 
is  similar  to  that  of  the  two-dimensional  theory  of  Cosserat  (or  directed) 
surfaces  which  are  relevant  to  fluid  sheets  and  to  shells.  For  a brief 
historical  account  of  the  developments  of  the  theory  of  directed  surfaces  and 
directed  curves  and  for  further  background  information  pertaining  to  a direct 
formulation  of  fluid  sheets  and  fluid  jets,  together  with  additional  references, 
the  reader  is  referred  to  a recent  expository  paper  by  Naghdi  [1].  The  first 
application  of  the  theory  of  a directed  curve  to  an  incompressible  Newtonian 
fluid  jet  was  given  by  Green  and  Laws  [2]  and  further  work  on  the  subject  was 
contained  in  a paper  by  Green  [3].  These  papers  are  concerned  with  a nonlinear 
theory  of  circular  jets,  which  includes  the  effects  of  both  surface  tension  and 
gravity.  More  recent  studies  on  the  subject,  which  will  be  referred  to  below, 
deal  with  temporal  instability  and  spatial  instability  of  incompressible  Newtonian 
viscous,  circular  jets  [4,5]. 

The  simplest  theory  appropriate  for  fluid  jets  that  can  be  constructed  on 
the  basis  of  a Cosserat  curve  comprises  a material  curve  and  a pair  of  directors. 
This  theory  is  suitable  for  many  applications  as  is  clearly  evident  from  the 
contents  of  papers  of  Caulk  and  Naghdi  [4]  and  of  Bogy  [5].  The  former  deals 
with  the  onset  of  breakup  of  a Newtonian  viscous  jet,  while  the  latter  is  mainly 
concerned  with  the  related  problem  of  jet  breakup  formulated  as  a boundary-value 
problem  in  connection  with  ink-jet  printing.  A brief  account  of  the  developments 
which  utilize  the  theory  of  a Cosserat  curve  with  two  directors  is  discussed  in 
section  3 of  this  paper. 
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3. 


Several  technologically  important  problems  in  rheology,  including  the 
swelling  effect,  are  concerned  with  fluid  flow  in  a pipe.  With  the  hope  that 
it  will  eventually  be  possible  to  treat  such  problems  by  direct  approach,  most 
of  the  remainder  of  the  paper  is  devoted  to  a Poiseuille  flow  which  requires  the 
use  of  the  next  hierarchical  theory  of  Cosserat  curves,  namely  that  comprising 
a material  curve  with  five  directors.  Instead  of  developing  separately  this 
next  hierarchical  theory  with  five  directors,  it  is  just  as  convenient  to  con- 
sider the  more  general  theory  of  Cosserat  curves  having  any  finite  number  of 
directors.  Thus,  in  section  2 we  construct  the  theory  of  directed  curves  with 
L (S2)  directors  and  then  briefly  discuss  the  results  of  the  special  theories 
when  the  number  of  directors  are  two  and  five,  respectively.  The  kinematics  of  a 
special  Poiseuille  flow  in  a straight  circular  pipe  is  considered  in  section  4 
in  the  context  of  an  approximation  procedure  whereby  the  position  vector  in  the 
three-dimensional  theory  is  taken  in  the  form  of  a Taylor  series  expansion  and  is 
then  assumed  to  be  quadratic  in  the  cross-section  coordinates  (see  Eq.  (Al8)  of 
Appendix  A).  The  kinematical  results  of  section  4 motivate  the  choice  of  the 
corresponding  kinematics]  ingredients  in  the  theory  of  a Cosserat  curve  with 
five  directors  The  latt  T is  employed  in  the  discussion  of  Poiseuill<  flow  in 
a circular  pipe  by  direct  approach  in  section  5.  It  should  be  noted  that  this 
solution  includes  the  ef !'■  ct  of  "normal  force"  corresponding  to  the  "normal  stress" 
effect  in  the  three-dimensional  theory.  In  fact,  the  relationship  between  the 
two  effects  is  evident  from  the  results  of  section  6,  where  the  identification 
of  certain  quantities  in  the  direct  theory  is  discussed. 


4. 


2.  Directed  fluid  .jets  with  L directors 

Deformable  media  which  are  modelled  by  a material  curve,  embedded  in  a 
Euclidean  3-space,  together  with  L (L§2)  directors  assigned  to  every  point  of 
the  curve  will  be  called  Cosserat  curves  or  directed  curves  and  may  be  con- 
veniently referred  to  as  R (K=l,2,...).  All  such  directed  curves  have  the 

same  material  curve  but  the  number  of  assigned  directors  differ  for  each^  ft  . 

K 

Let  c,  the  material  curve  of  in  the  present  configuration  at  time  t,  be 

defined  by  its  position  vector  r relative  to  a fixed  origin;  and  let  f be  a 

convected  (Lagrargian)  coordinate  identifying  points  along  the  curve.  Further, 

let  L directors  at  r be  denoted  by  the  vector  functions  d (a, , • • .a..  = 1,2 ; 

~ ~Qfla2’’*aN  1 2 TJ 

N = 1,2, . . . ,K) , which  are  assumed  to  be  symmetric  in  the  indices  Then, 


a motion  of  the  directed  curve  ft„  is  specified  by 


r = r(5,t)  , d _ = d 


~“lQ!2’”aN 


A 

d 

MalaZ' 


N 


(?,t)  , (N  = 1,2 , . . . ,K) 


The  velocity  and  the  director  velocities  are  defined  by 


v = r , w = d 

~ ~ ~Q!1«2--*q,n  ~Q!1q2'',“n 


(2.1) 


(2.2) 


where  a superposed  dot  denotes  material  time  differentiation  holding  5 fixed. 

Also,  the  tangent  vector  to  the  curve  c denoted  by  a„  is  given  by 

~o 

dr 

a0=a^(5,t)=— (|,t)  . (2.3) 

Corresponding  to  the  requirement  in  the  three-dimensional  theory  [see  (A7)  of 
Appendix  A)]  that  a nonzero  volume  cannot  be  continuously  deformed  into  a zero 
volume,  some  restrictions  on  the  directors  d are  necessary  but  we  leave 

~eW**°N 

these  unspecified  here.  It  is  more  convenient  to  specify  such  conditions  as 
they  arise  in  special  cases  of  the  general  theory. 

The  relationship  between  the  number  of  directors  L and  the  number  K in 
(2.1)  which  identifies  the  order  of  the  hierarchical  theory  of  Cosserat 


In  the  absence  of  the  directors,  we  merely  have  a one-dimensional  curve  which 
can  serve  as  a model  for  the  construction  of  string  theory  by  direct  approach. 
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curves  is  seen  to  be  L = s(N+l)  so  that 


|(Ktl)  + K 


(2.4) 


This  result  is  in  agreement  with  a corresponding  development  of  the  kine- 
matics of  a jet-like  (or  rod-like)  body  from  three-dimensional  theory  of  continuum 
mechanics  contained  in  the  paper  of  Green  et  al.  [6],  Green  and  Naghdi  [7]  and 
Green  et  al.  [8].  According  to  (2.4),  with  K=1  the  number  L = 2 and  the 
Cosserat  curve  R^ = R consists  of  a material  curve  and  a pair  of  directors 
attached  at  each  point  of  the  curve.  Similarly,  in  the  case  of  the  Cosserat 
curve  Rg,  the  number  of  directors  is  five  and  so  on. 

We  assume  that  the  kinetic  energy  per  unit  length  of  the  curve  c is  given  by 

K clcl  ...  Off. 

T = p [iv  • v + E y v * w 

~ ~ N=1  ~ "<YV“n 


K al’ ' ’aN®l" * '®M 


(2.5) 


+ t s y * w -w.  q J 

N=1,M=1  ~<VQiN  ^l”-eM 

where  p=  p(§,t)  is  the  mass  per  unit  length,  the  coefficients  y and 

• • • BivfP-i  • • • * al ' * ’ °N 

y are  functions  of  5 and  t,  y are  symmetric  with  respect  to 

* ’ * Pm  Pi  • ' * Pm®i * " " 

indices  c^.-.c^j,  y =y  and  are  also  symmetric  with  respect 

to  c^c^.  • -c^j  and  • "Pm’  In  the  sPecial  case  of  R1  with  L=2’  we  may  use  the  notations 


w = w 


a “lPl  aP 

y , y = yH 


(2.6) 


We  define  the  momentum  corresponding  to  the  velocity  v by 


ST  r ^ l ar**“N 

■r-  = p[v+  Z y v ] 

~ N=1  ~®1 * * * 


(2.7) 


per  unit  length  of  c.  Similarly,  momenta  corresponding  to  director  velocities 


ST  r V,0(N  . * “i’*-0(nPi,"Pm 

r--2 = p(y  v + 2 y . 

^ ...c^  M=1  ^1 


(2.8) 


per  unit  length  of  c. 


Although  the  coefficients  in  (2.5)  are  regarded  as  functions  of  both  5 and  t,  it 
will  be  proved  presently  that  they  are  in  fact  independent  of  t> 
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Consider  now  an  arbitrary  part  of  the  curve  c whose  element  of  arc  length  is 

X 

as  = (a  )2d5  , a = a • a . (2.9) 


Let  c be  bounded  by  and  %=%2  ( ^ < ^ ) and  define  the  following  quantities: 

The  contact  force  n = n(e,t)  and  the  contact  director  forces  p = 

r+s  /v#  * -< 


^1^2  * * * °*N 

P (£,t),  each  a three-dimensional  vector  field,  in  the  present  configura- 

tion^;  the  assigned  force  f = f(^,t)  per  unit  mass  and  the  assigned  director  forces 
ala2*'-aN  ala2'"aN 

jt  = l (§,t)  per  unit  mass;  the  intrinsic  (curve)  director  forces 

al“2 * ‘ * “N  alQ2*”aN,  s ala2  '“n  , 

n = n (|,t;  such  that  (a_„)  are  measured  per  unit  length 

~ ~ jj  ~ 

of  c and  make  no  contributions  to  the  balance  laws  for  moment  of  momentum  and 
energy;  the  specific  internal  energy  e = e(|,t);  the  specific  heat  supply  r = r(§,t) 
per  unit  time;  and  the  heat  flux  h=h(§,t)  along  c,  in  the  direction  of  increasing 
I,  per  unit  time.  The  contact  director  forces,  the  intrinsic  director  forces  and 
the  assigned  director  forces  are  all  assumed  to  be  symmetric  with  respect  to 


indices  cx-j^ ...  Q£^ • 


_2 

The  contact  force  n has  the  physical  dimension  of  [MT  ],  where  [M]  and  [T] 


stand  for  the  physical  dimensions  of  mass  and  time.  The  dimensions  of  the  vector 

“l°^-*-0(N 
fields  £ 


depend  on  the  choice  of  the  dimensions  for  the  directors. 


For  example,  if  each  of  the  vectors  d is  chosen  to  have  the 

~alQ2"  ,aM  0QU...0U  a1c^...aNP1P2...&M 

dimension  of  length,  then  the  coefficients  y , y 

0fj_®2  • • • 

in  (2.5)  are  dimensionless  and  the  vector  fields  p have  the  same 

dimensions  as  n.  A parallel  remark  applies  also  to  the  physical  dimension  of 

a2  * * ' °*N 

the  assigned  fields  l .In  the  present  paper  we  choose  the  directors 

to  have  the  dimension  of  length. 

In  terms  of  the  foregoing  definitions  of  the  various  field  quantities  and 
with  reference  to  the  present  configuration  at  time  t,  the  conservation  laws  of 


’Depending  on  the  physical  dimension  specified  for  the  directors,  the  field  quan- 
®1®2 • • 

titles  p may  also  be  referred  to  as  the  director  couples.  In  fact,  for 

M = 1 in  (2.1)  , corresponding  to  the  Cosserat  curve  ft  the  terminology  of  both 
d 1 ai  a 

the  director  force  vectors  and  the  director  couple  vectors  is  used  for  p = p 

a,o5'...afr 

in  the  literature.  A parallel  statement  holds  also  for  the  fields  l 
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the  purely  mechanical  theory  of  a directed,  curve  ft  are: 

K 


d rh  d r^2  V,#afo 

dt  J pds  = 0 ’ dt  ] py 

h h 

, K Qf,  • • • a.. 


ds  = 0 (N  = 1, . . . ,K)  , 


(2.10) 


|}2p[v+  Zy1  Nw  Ids  = J 2pfds  + [n]!2  , (2.11) 

dt  ~ N=1  ~°i*-'aN  ~ ~ h 


A l/L—*-**  * & 

^1  1 

d 1*^2  r * “i*”“n  l al”,0iN 

— pfrxv  + rxSy  w + E d x y v 

dt  J?1 N=1  ~®1  * " * aN  N=l~0fl""aN 

K a1...aN01...pM 

J'i:iv«xy  iLflRJds 

N=1,M=1  ~°1,'-Q(N  ~*V‘*M 


(2.12) 


[»5o  “l  " " “n 

= pfr  X f + d X l }ds 

' - ~ ~a1...oN  ~ 

®1  * * * ®M  So 

+ [r  xn  + d X P ].  > 

~ ~ ~°i-*'aN  ~ h 


where  use  is  made  of  the  notation 


(2.13) 


Cf]^  = f(52)-f(5i)  » 

The  first  two  of  the  above  equations,  namely  (2.10)1  and  (2.10)2,  represent  the 
mathematical  statement  of  conservation  of  mass  (or  inertia)  while  the  remaining 
equations  in  the  order  listed  are  mathematical  statements  of  the  conservation  of 
linear  momentum,  conservation  of  director  momenta  and  the  conservation  of  moment 
of  momentum  (including  contributions  from  both  the  ordinary  momentum  and  the 
director  momenta). 

We  also  record  here  the  law  of  conservation  of  energy  for  ft  , namely 
4r  [?2[pc  + T]ds  = [?2p(r  + f • v+  £ l 1 ^ . w )ds 

dt  J-  V ~ ~ N=l~  ~*L  — *S 


. r ai'*'“N  ,^2 

+ In  • v + E P • w -h] 

~ ~ N=l~  ~Qfl**-°N  h 


(2.1U) 


The  assigned  field  f in  (:!.10)  represents  the  combined  effect  of  (i)  an  integrated 
contribution  arising  from  the  three-dimensional  body  force  denoted  by  f , e.g., 

M) 


<-■-  - ■ , V.  - 


that  due  to  gravitational  acceleration,  and  (ii)  an  integrated  contribution  of 

the  stress  vector  on  the  lateral  surface  of  the  jet-like  body*  denoted  by  f . 

ala2  * * " °N 

A parallel  statement  holds  for  the  assigned  fields  l . Similarly,  the 

assigned  heat  supply  r represents  the  combined  effect  of  (i)  an  integrated 
contribution  arising  from  the  three-dimensional  heat  supply  denoted  by  r^,  and 
(ii)  an  integrated  contribution  of  the  heat  supply  entering  the  lateral  surface 
of  the  jet-like  body  from  the  surrounding  environment,  denoted  by  r^.  Thus,  we 
may  write 


f = 


f.  +f 
~b  ~c 


“l“2-‘*%  “la2-*-V 

£ = £b  +*c 


r = r,  + r 
b c 


(2.15) 


The  various  quantities  in  (2.15)  are  free  to  be  specified  in  a manner  which 

depends  on  the  particular  application  in  mind  and,  in  the  context  of  the 

theory  of  Cosserat  curves,  the  inertia  coefficients  in  (2.4)  and  the  mass 

density  p require  constitutive  equations.  Indeed,  f , l and  r 

Q(l®2  * * * 

(or  certain  of  their  features),  as  well  as  and  r^,  can  be 

identified  with  the  corresponding  expressions  in  a derivation  from  the  three- 
dimensional  equations.  Likewise,  the  inertia  coefficients  in  (2.4)  and  the 
mass  density  p may  be  identified  with  easily  accessible  results  from  the 
three-dimensional  equations. 

Assuming  smoothness,  we  may  deduce  the  local  form  of  the  balance  equations 

(2.!°)i,2  to  be 


X - 0 


(N  = 1 , . . . ,K) 


1 

X = X(?)  = pa^3  , 


(2.16) 


*A  definition  for  jet-like  bodies  is  provided  in  Appendix  A.  Its  lateral 
surface  is  specified  by  equation  (Al4). 


are  independent  of  t and  are  functions  of  5 only.  Then, 


®1 • • • 

so  that  X.  and  y 

the  remaining  conservation  laws  (2.11)  to  (2.13)  yield  the  field  equations 

dn  _ 

7-r  + \f  = ° , 

05  ~ ~ 

* al * * " °N 

SP  al'--aN  ai * * * °N 

+ Xq.  — tt  > (N— j 


a_  x r -r  Z ( d 


N=1  - - * aN 


dd 

ai ’ ' ' aN  . ~aVaN  “i,,,an1 
x n + ^ x p > - 0 


a? 


(2.17) 

(2.1 8) 

(2.19) 


where 


K . 

f = f - v - z y w 9 (n  = 1 , . . . ,k)  , 

~ ~ ~ N=1  ~®1 ' * ' °N 


al‘"-°N  al * * ' aN • K al’**aN^l — • 

1 = i -y  v-zy  o 

- ~ ~ M=1  ~Pi#”Pm 


(2.20) 

, (N  = 1, . . . ,K)  . (2.21) 


With  the  help  of  (2.l6)  to  (2.19),  the  equation  for  conservation  of  energy  can 
be  reduced  to 


dv  K a-,* 
P = n • -r+  Z tt 
~ N=1 


Xr-Xe-fJ+P  = 0 , 

os 


,,  dw 

K V-On  ~®i’  ’ 
w + £ p rr- 

n=i~ 


(2.22) 

(2.23) 


where  P is  the  mechanical  power.  It  may  be  noted  here  that  within  the  scope  of 
the  purely  mechanical  theory,  the  expression  (2.23)  for  mechanical  work  can  be 
obtained  also  by  considering  the  rate  of  work  by  all  contact  force  and  director 
forces  and  all  assigned  force  and  assigned  director  forces  acting  on  the  curve  c 
and  its  end  points  minus  the  rate  of  increase  of  the  kinetic  energy  and  by  setting 

■»§o 

this  equal  to  | P d|;  in  this  connection,  see  for  example  [1,  Eq.  (3-l6)]. 

JSl 

In  the  remainder  of  the  paper  we  shall  be  concerned  with  special  cases  of 
the  above  equations  appropriate  for  directed  curves  ft-  and-  1°  particular, 


for  the  directed  curve  9.^,  the  appropriate  differential  equations  of  the 


! ' * 
i ! 


mechanical  theory  consist  of  (2.16)2  with  N-2,  (2.16)^  (2.17)  and. 


..  a c* 

+ \q  = tt  , 


Zr  + n+.J*  , 

a xn  + d x tt“  + d x tt°^  + X p“  + -rf®’  X P°^  = 0 
~3  ~ ~ ~ap  ~ 05  ~ ~ 


Also,  the  expression  for  the  mechanical  power  in  this  case  is  given  by 

dv  „ 3w  3w 

P ..  if  • 

r-**  35  <^-0(  ^ ~ap  ^ ^ o 5 


(2.24) 

(2.25) 

(2.26) 

(2.27) 


A special  case  of  the  above  general  development  with  two  directors,  i.e., 
for  a Cosserat  curve  R=R1}  was  first  given  by  Green  and  Laws  [9]  in  the  context 
of  thermomechanics.  A related  development  of  a mechanical  theory  employing 
three  directors  is  contained  in  a paper  by  Cohen  [10].  Further,  aspects  of  the 
basic  theory  with  two  directors  appropriate  for  a Cosserat  curve  R are  contained 
in  a more  recent  paper  by  Green  et  al.  [11]. 


3.  Straight  Newtonian  viscous  .jets. 

We  consider  here  a special  case  of  the  general  theory  of  section  2 for  a 
directed  curve  = ft , which  may  include  the  effects  of  gravity  and  surface 
tension,  and  go  on  to  illustrate  its  applicability  to  certain  Newtonian  flows 
and  related  problems  of  instability  or  breakup  of  straight,  incompressible, 
viscous  jets  of  circular  cross-section.  The  motion  of  the  directed  curve  ft  is 
specified  by  (2.1)^  and  by  (2.1)^  with  N=l.  Also,  the  velocity,  the  director 
velocity  and  the  inertia  coefficients  in  this  case  are  given  by  (2.2)^,  (2.6)^ 

and  (2.6)  respectively.  The  only  kinetical  quantities  which  occur  in  the 

* & & 

theory  of  a directed  curve  with  only  two  directors  are  the  forces  n,p  ,rr  and 


the  fields  f,i  . Thus,  in  the  context  of  the  purely  mechanical  theory,  the 
local  forms  of  the  relevant  conservation  laws  are  the  mass  conservation 
(2.16)  , the  equations  of  motion  (2.17)  and  (2.24)  and  the  consequence  of 

1 JO 

moment  of  momentum  is  given  by  (2.26)  after  omitting  terms  which  involve  d 

~ap 

and  3d  /d£. 

~c*P 

For  a constrained  theory  of  the  directed  curve  ft,  we  assume  that  each  of 
the  functions  n,rra,pa  is  determined  to  within  an  additive  constraint  response 
so  that 

_ a rof.Aa  -a  , 

n = n + n , n = tt  + tt  , p = p+p  , 

A Ajy 

where  n,rr  ,£  are  determined  by  constitutive  equations  and  the  constraint 
responses  n,^,^1  are  arbitrary  functions  of  £,t  and  do  no  work.  For  the  class 
of  fluid  jet  problems  discussed  in  [3>4,5],  the  condition  of  incompressibility 
can  be  shown  to  yield 


(3-D 


dt  ° 


(3.2) 


/ 


* 


Then,  assuming  that  for  an  incompressible  Newtonian  viscous  fluid  at  constant 


These  are  special  cases  of  those  defined  in  section  2 following  (2.9). 


12. 


• ...  . ■ . >.  •■*  '•»  J.  v.f  '-n  . . : 
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temperature  the  constraint  response  n,TT  ,p  do  not  depend  explicitly  on  the 


kinematical  quantities 


dv  dw 

~ ~<v 

— - T*r  — « 

d?  * lia  * S5  ’ 


(3.3) 


it  can  be  shown  that  [3,12] 


n = - qd  x dp  , rr“  = - qe“Pd^  X a , pa  = 0 , 


(3.4) 


where  the  Lagrange  multiplier  q is  an  arbitrary  scalar  function  of  |,t  and  ec 
is  defined  by 


11  22  _ 12  21 
e =e  = 0 , e = - e = 1 


(3.5) 


We  now  specialize  the  results  of  the  previous  section  to  straight  jets 
of  elliptical  cross-section.  In  order  to  display  some  details  of  the  kine- 
matics of  a straight  jet,  including  the  rotation  of  the  directors  in  a plane 
normal  to  the  jet  axis,  it  is  convenient  to  introduce  a fixed  system  of 
rectangular  Cartesian  coordinates  (x,y,z)  with  the  z-axis  parallel  to  the 
jet.  Further,  let  the  unit  base  vectors  of  the  rectangular  Cartesian  axes 
be  denoted  by  (i,j,k)  and  introduce,  for  later  convenience,  the  additional 
base  vectors 


e = i cos  e + j sin  0 , e0 


= -i  sin  0 + j cos  0 , e_  = k , 


(3.6) 


where  0 is  a smooth  function  of  z and  t.  We  assume  that  the  directors  are  so 

restricted  that  they  describe  an  elliptical  cross-section  of  smoothly  varying 

orientation  along  the  length  of  the  jet  and  that  at  each  z= const.,  the  base 

vectors  e^  and  e^  lie  along  the  major  and  minor  axes  of  the  ellipse,  recpec- 

A, 

tively.  Then,  the  angle  0=0(z,t),  called  the  sectional  orientation,  specifies 
the  orientation  of  the  cross-section  at  z = constant  as  a function  of  time.  With 
this  background,  we  new  restrict  motions  of  the  directed  curve  R such  that  in 


,{ 


i 


the  present  configuration  at  time  t, 


£ = £(5»t)  , da_  * 0!^  , dg  = 02eg  , 

01  = £].(?, t)  , 02  = 02(?,t)  , 

where  0^  and  02  measure  the  semiaxes  of  the  elliptical  cross-section.  In  the 
case  of  a circular  jet,  0^  = 02* 

The  complete  theory  of  the  directed  jet  of  this  section  also  requires  the 

specification  of  explicit  values  for  X,ya,y°^,f  and  ia.  In  particular,  the 

values  for  X,ya,y® ^ may  be  obtained  by  an  appeal  to  certain  results  from  the 

three-dimensional  description  of  the  jet.  For  this  purpose,  we  may  approximate 

Ax  . 

the  vector  function  r on  the  right-hand  side  of  (Al)  of  Appendix  A and  write 

rw 

t. 

this  as 


(3-7) 


r = r + 0 d 

Further,  we  choose  the  curve  0® = 0 as  the  line  of  centroids  of  the  jet-like 
body  and  identify  this  curve  with  the  curve  c in  the  theory  of  a Cosserat 
curve.  This  leads  to  the  identification 


(3.8) 


X = P(a33)2  = J^P  g2d01 


de2  , 


xy' 


M- 

Ja 


P g20Qrd0xde<"  = 0 , Xy®0  = f p*g20®0Kd0~d0“ 

•»  r* 


J , „1J  2 


rr2 


(3.9) 


where  p is  the  three-dimensional  mass  density  in  (All)  and  the  determinant  g 

defined  by  (A5)^  is  calculated  from  the  approximation  (3*8) . Again,  with  the 

use  of  (3*7)  and  the*  equations  of  motion  (2.9)0  the  expressions  for  f and  la 

* 

can  be  identified  in  terms  of  the  integrated  body  force  f over  the  cross- 
section  a and  specified  pressure  and  surface  tension  over  the  boundary 
3d  of  d (for  details,  see,  for  example,  Caulk  and  Naghdi  [12]).  We  observe 
that  since  y®  = 0 by  (3.9)2>  the  equations  of  motion  (3.8)  and  (3-9)  assume  a 
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slightly  simpler  form.  We  do  not  record  here  the  system  of  ordinary 
differential  equations  which  can  he  obtained  for  both  inviscid  and  linear 
viscous  fluids.  These  are  readily  available  elsewhere  and  are  also  dis- 
cussed  in  [1], 

The  development  of  constitutive  equations  for  an  incompressible  Newtonian 
viscous  jet  of  both  circular  [3]  and  elliptical  [4]  cross-sections  are  readily 
available  in  the  papers  cited,  and  the  constitutive  coefficients  are  identified 
in  terms  of  the  shear  viscosity  and  the  geometry  of  the  cross-section  of  the 
jet.  We  do  not  elaborate  here  on  these  results  and  refer  the  reader  to  the 
papers  of  Green  C 3 J and  Caulk  and  Naghdi  [4]. 

In  the  rest  of  this  section  we  briefly  describe  some  evidence  of  the 
relevance  and  applicability  of  the  direct  formulation  of  viscous  fluid  jets, 
especially  to  problems  of  instability  of  viscous  jets  which  utilize  the  basic 
equations  of  the  theory  of  direct  fluid  jets  of  this  section.  For  definite- 
ness we  consider  the  linearized  version  of  the  relevant  equations,  neglect  the 
effect  of  gravity  and  discuss  the  onset  of  breakup  of  a viscous  jet  due  to 
surface  tension,  i.e.,  the  so-called  capillary  instability.  Although  our 
interest  here  centers  mainly  on  Newtonian  viscous  jets,  in  order  to  assess 
the  nature  of  the  prediction  of  the  direct  approach,  it  is  desirable  to  con- 
sider also  the  breakup  of  an  inviscid  jet  since  this  enables  us  to  compare 
the  results  with  the  available  exact  analysis  of  the  breakup  of  inviscid  jets 

* 

due  to  Rayleigh  obtained  by  means  of  the  linearized  three-dimensional  equations. 
For  both  inviscid  and  viscous  jets,  Rayleigh  derived  the  explicit  result  that 

Ithe  jet  is  unstable  only  in  the  axisymmetric  mode  of  disturbance.  Inasmuch 

as  the  direct  theory  considered  here  does  not  begin  with  the  three-dimensional 
equations,  all  modes  of  disturbance  which  occur  in  the  present  one -dimensional 

# 

References  to  Rayleigh's  papers  on  the  subject  are  cited  in  [4]. 
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direct  theory  must  be  examined  for  stability.  Thus,  in  order  to  allow  for 
the  growth  of  a general  disturbance  which  is  not  necessarily  symmetric, 

Caulk  and  Naghdi  utilizing  the  results  of  [12]  derive  in  [4]  a system  of 
linearized  equations  governing  the  small  motions  of  a (nonrotating)  incom- 
pressible inviscid  jet  of  elliptical  cross-section  superposed  on  uniform 
flow  of  a circular  jet.  They  show  that  the  solution  to  these  linearized 
equations  can  be  decomposed  into  two  modes,  representing  a symmetric  and  an 
anti -symmetric  disturbance  in  the  shape  of  the  free  surface.  The  anti- 
symmetric mode  is  stable  for  all  wavelengths,  while  the  symmetric  mode  is 
found  to  be  unstable  over  a range  of  longer  wavelengths.  In  terms  of  a 
description  of  growth  in  the  unstable  mode,  comparison  of  the  conclusions 
is  found  to  agree  extremely  well  with  the  corresponding  exact  three- 
dimensional  analysis  of  Rayleigh. 

In  the  case  of  a straight  incompressible  Newtonian  viscous  jet,  through 
a comparison  with  available  three-dimensional  numerical  results  (Chandrasekhar 
[13]),  the  solution  obtained  is  shown  to  be  an  improvement  over  an  existing 
approximate  solution  of  the  problem  by  Weber  [l4].  A related  study  by  Bogy 
[5],  concerning  the  instability  of  an  incompressible  viscous  liquid  jet  of 
circular  section,  partly  overlaps  with  the  work  of  Caulk  and  Naghdi  [4]  on 
the  temporal  instability  of  a viscous  jet.  Bogy  [5]  confines  attention  to 
the  symmetric  mode  of  disturbance,  and  considers  mainly  the  spatial  instability 
of  a semi-infinite  jet  formulated  as  a boundary-value  problem.  For  additional 
background  on  breakup  and  drop  formation  in  viscous  jets  of  circular  cross- 
section,  see  a recent  article  by  Bogy  [15]  which  contains  additional  references 
on  the  subject. 

* 

If,  in  the  context  of  the  direct  theory  of  this  section,  the  stability  analysis 
is  confined  through  £ priori  assumptions  to  the  symmetric  mode  of  disturbance 
only,  then  any  conditions  for  instability  are  only  sufficient. 


% 
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4 . Kinematics  of  flow  in  a pipe:  A special  Poiseuille  flow. 

We  elaborate  here  on  some  three-dimensional  kinematical  aspects  of  a fluid 

flow  in  a straight  circular  pipe,  using  the  notation  of  Appendix  A.  We  first 

recall  that  for  a Poiseuille  flow  in  a fixed,  infinite,  circular  pipe,  the  only 

* 

nonvanishing  component  of  the  velocity  vector  v is  its  axial  component  and  this 
is  a function  of  the  radius  of  the  cross-section.  Moreover,  the  path  lines  are 
straight  and  parallel  to  the  axis  of  the  pipe;  and  the  boundary  conditions 
imposed  on  the  flow  are  the  vanishing  of  the  axial  velocity  at  the  wall  of  the 
pipe  and  continuity  of  the  stress  tensor  at  the  center  line  of  the  pipe.  We 
keep  this  background  information  in  mind  in  the  development  of  this  section. 

Consider  now  the  kinematics  of  flow  in  a pipe  of  uniform  circular  cross- 
section  in  the  context  of  an  approximation  procedure  in  which  the  position 
vector  is  approximated  by  the  expression  (Al8)  of  Appendix  A.  Let  x^ = (x,y,z) 
be  a fixed  right-handed  rectangular  Cartesian  coordinate  system  and  let  the 
associated  unit  base  vectors  be  denoted  by  i,j,k  (see  Fig.  2).  We  identify  the 
center  line  of  the  pipe  with  the  z-axis  and  specify  its  cross-section  of  constant 
radius  a by 


2 2 2 
x +y  = a 


(4.1) 


It  is  convenient  to  consider  first  the  kinematics  of  the  flow  at  time  T§t. 

*.  . A*,  i . . A 

Thus,  recalling  the  rotation  of  Appendix  A,  let  r (r)=r  (0  ,t),  r (t) = r (| , t) , 

d (-r)  = d (§,t),  d (Tj  = d (£,t)  designate  the  various  quantities  in  (Al8)  at 
'-op  '-orp 

time  t;  and,  with  reference  to  the  present  configuration  at  time  t,  we  adopt 
the  notations  r =r  (t),  r = r(t),  d =d  (t),  d =d  _(t).  Then,  given  the 

***  r>-*  rs9  ~-c*p 
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approximation  (Al8)  the  position  vector  r (t)  and  the  velocity  v (t)  are: 


(t)  = r(t)  + 9°d^(T)  + , 


(4.2) 


(t)  = v(t)  + 9°W  (t)  + 909^w  (t)  , 


(4.3) 


where  v(t),  w (t),  w _(t)  are  defined  as  in  (A19).  The  position  vector  of  the 
~ ~a  ~ap 

center  line  of  the  pipe  and  the  velocity  of  a particle  on  the  center  line  are, 
respectively,  r (0,0, e,r) = r(f ,t)  and  v (0,0, ?,t) = v(t) . Since  the  z-axis  is 
taken  to  he  coincident  with  the  center  line  and  since  in  a Poiseuille  flow 
particles  on  the  center  line  move  with  a constant  velocity,  we  set 


t(t)  = z(i)k  , v(t)  = z(-r)k  = vk  , 

rsj  rsj  rsj  rw 


(4.4) 


where  z(t) = z(0,0,|,t),  a superposed  dot  in  (4.4)  and  elsewhere  in  this  section 

designates  material  time  differentiation  with  respect  to  t and  v is  a positive 

constant.  At  this  point,  without  loss  in  generality,  we  identify  the  convected 

coordinate  f of  a particle  with  its  z-coordinate  in  the  reference  configuration 

at  time  t , i. e. , 
o ’ 


l = *(t0)  = z(e1,e2,?,T0) 


(4.5) 


Hence,  for  particles  on  the  center  line  of  the  pipe  in- the  reference  configura- 
tion we  have 


§ = z(tq)  = 2(0,0, 5, Tq)  , 

r(T  ) = z(t  )k  = §k  . 

/■»«>  o O ~ 


(4.6) 
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We  also  choose  the  convected  coordinates  0 and  0 such  that  they  are,  respec- 
tively, — times  the  x and  y coordinates  of  a particle  in  the  reference  configuration 
8 

at  time  t . Moreover,  observing  that  the  x,y  coordinates  of  a particle  do  not 
change  in  the  special  flow  under  consideration,  we  set 


91  = a°  = > x(t)  = xf©1,©2,^)  » 

, y(T)  - yA(01,02,?,T)  . 


(4.7) 


Integration  of  the  second  equation  in  (4.4)2  and  the  use  of  (4.6)^  results  in 


' 


z(t)  = v(t-to)  + z(tq)  = v(t-to)  + ? , 


(4.8) 


so  that  (4.4)^  may  now  be  written  at 


r(-r)  = r (t  ) + v(t-t  )k 

O O r* 


= [?  + v(T-TD)]k 


(4.9) 


In  order  to  obtain  some  restrictions  on  the  functions  d (t)  and  d (t) 

~ap 

in  the  expression  for  r (t),  consider  first  the  scalar  product  of  (4.2)  and 

rw 

1 • G • , 

t*(t)  • k = z(t)  + 01d1 (t)  • k + 02d  (t)  • k 

+ (e1)2d1i(T)  • k + (e2)2d22(T)  • k + 26192d12(T)  • k . (4.10) 

* 

In  view  of  the  symmetry  of  the  flow,  the  scalar  r • k must  remain  unaltered 
under  the  transformations 


(a)  x--y,  y-y  ; (b)  x-x,  y-  -y  ; (c)  x--x,  y-  -y  . 


(4.11) 


Hence,  from  the  symmetry  transformations  (4.11),  in  the  order  listed,  we  obtain 

e1d1(r)  • k + 20102d  • k = 0 , 


92d2(T)  • k + 20102d12  • k = 0 , 


(4.12) 


0^1,  (t)  • k + e2d_(T)  • k = 0 . 


It  follows  from  the  above  three  conditions  that 


<1,5(1)  • k = 0 , d (t)  • k = 0 , dp(T)  • k = 0 . 


(4.13) 
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According  to  (4.13) , the  vector  functions  d (t)  and  d10(r)  have  no  components 

along  the  axis  of  the  pipe  and  thus  cannot  contribute  to  the  axial  velocity 

-#■ 

v (t).  Thus,  in  view  of  (4.13)0  n>  we  set 
~ >5 

d (t)  = ai  , <U(t)  = aj  . (4.14) 

rwcL 

Next,  using  (4.4).  end  (4.l4),  from  the  scalar  product  of  (4.2)  and  i we  have 

_L  rsj 

r*(T)  • i - 91a  + (e1)^  (t)  • i + (02)V9(t)  • i + ae^-e^  • i . (4.15) 

* 

In  view  of  the  symmetry  of  the  flow,  the  scalar  r • i must  remain  unaltered 

/V  A/ 

under  the  transformation  (b)  in  (4.1l)  and  the  resulting  expression  when 

1 2 

compared  with  (4.15)  yields  40  0 d.n  • i = 0 so  that 

r*d.d 

^(t)  • i = 0 . (4.16) 

* . . 

Similarly,  by  requiring  that  r • j be  unaltered  under  the  transformation  (a) 

Arf  /v 

in  (4.11),  we  arrive  at 

~12 ( T ) ’ A = 0 • (4,17) 

It  follows  from  (4.13)^ 5 (4.16)  and  (4.17)  that 

d.„(T)  = 0 . (4.18) 

~Ltr.  r*0 

On  substitution  of  (4.9) > (4.14)  and  (4.18),  (4.2)  and  (4.3)  reduce  to 

r*(t)  = v(T-T0)^,+  r(T0)  + a(01i+  02^)  + (e1)2^.^)  + (S^dggCf)  (4.19) 

and 

v#(t)  = vk+ (e1)2w11  + (e2)2«g2  , (4.20) 

Again  the  symmetry  of  the  flow  requires  that  at  points  along  the  x-axis 
and  the  y-axis  we  have  v (0,1,5,t)  = v (1,0, 5,t).  This  leads  us  to  conclude  that 


■ 

_ 


20. 


' ^2(t) 


(4.21) 


and  (4.20)  can  be  rewritten  as 


v (t)  = vk+ [(01)2  + (02)2jw^;l(t)  . 


(4.22) 


From  the  condition  that  the  particle  velocity  vanishes  at  the  wall  of  the  pipe, 


i.e. , 


v*(t)  = 0 when  (01)2  + (02)2  = 1 , 


(4.23) 


we  have 


w (t)  = - vk  . 
'xLJ.  *** 


Integration  of  (4.23)  with  respect  to  t yields 


(4.24) 


^(t)  = - v(T-TQ)k  , (4.25) 

where  the  reference  value  of  ^-q(t)  at  tq,  i.e.,  <I^(to),  arising  from  the 
integration  has  been  set  equal  to  zero.  Recalling  (4.20),  in  a similar  manner 


we  also  obtain 


2P2(T)  = £l1(t)  =-v(T-T0^  ‘ 


X22  XU"'  o'~ 

We  now  summarize  the  foregoing  results  and  record  below  the  expressions 
for  the  position  vector  and  the  particle  velocity  both  in  the  reference 
configuration  at  time  tq  and  at  time  t: 

r *(t  ) = ?k  + a(01i  + 02j)  = §k  + xi  + yj  , 


(4.26) 


(4.27) 
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r (t)  - a^i  + 02j)  + {§  + v(t-t  )[1  - (Q1)2  - (02)2]}k 


“ xi,+yJ  + (?  + v(T-To)[1‘^2_^]5k 


a a 


' re  +{5  + v(t-t  , 


(4.28) 


f(r)  = [i  - (e1)2  - (e2)2]vk  = (1  - (f)2  - 02]vk  = [1  - 02]vk  , 


(4.29) 


2 2 - 

where  r=  Tx  +y  Y is  the  radial  distance  in  cylindrical  polar  coordinates  (r,8,z) 
with  associated  orthonormal  base  vectors  (e  ,e  ,k).  It  is  clear  from  (4.28)  that 


:(t)  = § + v(t-t  )(1-^)  , 


(4.30) 


which  reduces  to  (4.8)  on  the  center  line  of  the  pipe  (r  = 0).  It  is  clear  from 

y 

(4.29)^  that  (dv  (t)/St)  = 0 for  fixed  values  of  x and  y and  hence  the  motion  is 

y. 

steady.  Also,  using  (4.29) 5 we  may  calculate  the  rate  of  shear  K as  ordinarily 
defined  for  helical  flow  in  the  context  of  the  three-dimensional  theory  (see,  for 
example,  [16,  p.  23]).  Thus,  since  the  expression  (4.29)^  i6  a function  of  r 
alone,  we  have 


dv*(T) 

I-&H 


2v 


(4.31) 


where  v_(t)  = v(t)  • ic  and  the  double  vertical  bar  stands  for  the  absolute  value. 
3 ~ 

* y- 

Figures  2 and  3 depict  the  position  vectors  r (t)  and  r (t  ) for  the 
special  Poiseuille  flow  in  a circular  pipe  discussed  in  this  section.  These 
figures  explicitly  indicate  the  interpretations  associated  with  the  functions 

d (t)  and  d „(t)  as  given  by  (4.l4),  (4.l8)  and  (4.26).  It  can  be  seen  from 

~o  ~ap 

Figs.  2(b)  and  2(c)  that  fluid,  which  in  a reference  configuration  at  time  t0 
occupies  a cylindrical  region  such  as  that  bounded  by  the  normal  sections  at 
5 and  l,  in  the  configuration  at  time  r occupies  the  region  between  the 
paraboloids  of  revolution  CA 'D  and  EB'F.  Figure  3 exhibits  how  a material  line 
such  as  BE  in  the  configuration  at  time  deforms  into  a portion  of  a parabola 
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indicated  by  B'E. 


Before  closing  this  section,  we  indicate  that  the  flow  characterized  by 
(4.28)  is  isochoric,  i.e.,  the  volume  of  each  part  of  the  body  remains  unchanged. 

With  the  use  of  (4.28)^,  (A5)1  and  the  notation  (A2)^  of  Appendix  A,  the 

— * 
scalar  g2  which  occurs  in  (A7)  is 


g2  = & * fc  * g3 


which  is  a constant. 


Hence 


_d_ 

dT 


(4.32) 


T 


k- 


5.  A Poiseuille  flow  in  the  theory  of  directed  jets. 

We  discuss  here  a Poiseuille  flow  for  an  incompressible  viscoelastic  fluid 
in  a circular  pipe  employing  the  theory  of  a directed  curve  R^  in  the  absence 
of  the  effects  of  gravity  and  surface  tension.  The  motion  of  R^  is  specified 
by  (2.1)^  and  (2.1)2  with  K=2.  It  is  convenient  to  consider  first  the  motion 
of  a't  time  t 35 1.  Thus,  we  write 


r(T)  - f(?,T)  , ^(t)  = ^(5,t)  , ^(t)  = ^p(?5T) 


(5.1) 


and,  with  reference  to  the  present  configuration  at  time  t,  we  adopt  the 
notations 


r = r(t)  , d = d (t)  , h.  = d (t)  . (5-2) 

~ ~ ~a  ~a  ~a(3  ~ap 

As  in  section  4,  again  let  x^  = (x,y,z)  be  a fixed  right-handed  rectangular 
Cartesian  coordinate  system  but  here  we  use  the  notation  e^  for  the  unit  base 
vectors  in  place  of  i,j,k. 

rsj 

Let  the  curve  c of  R^  be  the  center  line  of  the  pipe,  and  identify  the 
latter  with  the  z-axis.  Also,  we  choose  the  directors  d (t)  such  that  they 
describe  the  symmetries  of  the  flow.  In  a Poiseuille  flow  the  path  lines  are 
straight  and  parallel  to  the  axis  of  the  pipe,  the  velocity  of  the  center  line 
is  constant  and  the  particles  located  at  the  wall  of  the  pipe  must  remain 
there.  With  this  brckground  and  guided  by  the  kinematical  results  of  section  4, 
for  the  motion  under  consideration  we  set 


t(t)  = r(r  )+v(t-t  )e  , r(r)  = , d (t)  = ae_  , 

~ ~ O O -O  — O " 


-3  » V 'o'  *~3 

djii(T)  = ^(t)  =-v(T-To)e3  , ^ = 0 , 


(5-3) 


where  r(TQ)  is  the  position  vector  of  the  fluid  particle  on  the  center  line  in 
some  reference  configuration  at  time  tq  and  it  should  be  noted  that  for  the 
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motion  (5*3)  the  tangent  vector  a = e_.  By  (2.2),  the  velocity  v and  the 
director  velocities  at  time  t are 


I(t)  = v%  » ISll  = ^2  =_v£3  5 


(5-4) 


^ = & = *02  = ° • 


With  the  above  choice  of  kinematical  ingredients,  it  is  clear  that  the  motion 
is  steady. 

To  complete  the  theory  of  the  directed  curve  for  an  incompressible  medium 
(see  the  discussion  at  the  end  of  Appendix  A following* (A26) ) , we  assume  that 


each  of  the  functions 


a a cyB  a& 
n , n , p , tt  p 


in  (2.27)  is  determined  to  within  an  additive  constraint  response  so  that 

a -a.'V  a ra  , (Nx 
n = n + n , n = tt  + rr  , p = P+P  , 

r+j  /v  ^ (V  rw  r+s 


(5.5) 


(5.6) 


where 


A /V  /V  Aap  /Vp 

? ! ) P 5 P 


(5-7) 


are  determined  by  constitutive  equations  and  the  constraint  responses 


- -a  ra  -orP  rap 
n , tt  , p , rr  , P 


(5.8) 


are  arbitrary  functions  of  §,t  and  do  no  work.  We  assume  that  the  constraint 

responses  (5*8)  do  not  depend  explicitly  on  the  gradient  of  the  velocities 

v(t),  w (t),  w (t)  and  since  they  are  workless  we  have 
~ ~a  ~orp 
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(5-9) 


l ! 


n • v ' + ir5*  • w +naP.w  +ia.w'  + FaP-w'  = 0 
~ ~ ~cr  ~aP  ~Qf  ~afi 

where  for  convenience  we  have  introduced  the  notation 

< )'.£<>  • 

Next,  we  introduce  a set  of  Lagrange  multipliers 

q , qY  , qY71  , ^ 


(5.10) 


(5.11) 


and  note  that  each  of  the  multipliers  q,  = (q\q^),  q^=  (q^,q^,q^,q^) , 
etc.,  may  depend  on  5 and  t.  Then,  incorporating  the  conditions  (A28)  to  (A3l) 
of  Appendix  A in  (5*9)  by  the  usual  procedure,  after  a lengthy  manipulation  we 
obtain  the  appropriate  expressions  for  the  constraint  responses  (5*8)  in  terms 
of  (5.U  ) and  the  kinematic  quantities  (5.2)  and  their  spatial  derivatives. 

Since  some  of  the  expressions  for  the  constraint  responses  are  quite  lengthy 
and,  in  any  case,  will  not  be  needed  in  the  present  development,  we  do  not 
record  them  here.  We  note,  however,  a typical  expression,  namely 

n =-e^q^x^  + q\x^p+2qy\llBxdYp]  , q^  = q^  . (5.12) 

In  view  of  the  choice  (5*3)  > most  of  the  terms  in  the  constraint  responses 
vanish  for  the  special  Foiseuille  flow  under  consideration;  and,  in  particular . 
the  expression  (5*12)  reduces  to 

n = - qa2^  > av(T-To)qVe^  . (5-13) 

The  response  of  an  incompressible  viscoelastic  directed  fluid  jet  may 

A A A , 

depend  on  the  entire  history  of  the  functions  r,d  ,d  _ in  (5-1)  and  their 
* 

derivatives.  As  in  the  corresponding  constrained  three-dimensional  theory,  we 
assume  that  the  quantities  (5*7)  are  determined  by  the  history  of 

U(T)  = (a_(T),d  (T),dfl(T),d/(T),d'ft(T))  . 

~or  ~ap  ~a  ~Cfp 


(5-14) 


AAA 


While  the  derivatives  may  be  calculated  from  the  histories  of  r,c^,d. 
if  the  latter  are  sufficiently  smooth,  their  inclusion  is  for  explicitness. 
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that  for  the  flow  under  discussion 


It  is  clear  from  (2.3)  and  (5-3)0  , . 

2 ,3>4 

the  above  kinematical  quantities  reduce  to 


b(-r)  = (e  ,ae  ,-v(t-t  )e  ,0,0)  . 

~o  ~ct  o ~3 

A 

A typical  constitutive  equation,  for  example  that  for  n,  may  be  written  as 


(5.15) 


n=£(v(t))  . (5-16) 

-00 

But,  since  the  argument  v(t)  of  the  functional  in  (5.16)  has  the  form  (5.15), 
it  follows  that  the  right-hand  side  of  (5.16)  can  be  regarded  as  a function  of 
v and  we  have 

A A 

(5.17) 


n = n(v) 


with  similar  expressions  for  all  other  kinematical  quantities.  We  require  that 
all  constitutive  functionals  or  constitutive  functions,  such  as  those  in  (5.16) 
and  (5.I7),  vanish  for  v = 0.  This  requirement  is  similar  to  the  normalization 
of  the  constitutive  response  functionals  in  the  three-dimensional  theory. 

Referred  to  the  orthonormal  basis  e.  = (i,j,k),  the  functions  (5.5)  can  be 

rsJ  rs*  r*f 

expressed  in  terms  of  their  components  in  the  form 


Jo  - “i£i 

ap  _ op 


= tr 


i~i  » & " •'•iSi 

oP  aP 
i~i  ~ * • i~i 


(5.17) 


Also,  since  the  effect  of  gravity  is  neglected,  f^  = 0 , £^=0,  A^=0  anc^  ^he 


notation  of  (2.15),  we  have 

1 yd 


f = f , e - «“  , = g?  . 

The  quantities  (5.18)  may  be  regarded  as  representing  the  forces  that  must  be 
supplied  by  the  pipe  wall  in  order  to  sustain  the  assumed  flow. 

As  a consequence  of  the  symmetry  of  the  assumed  flow,  an  examination  of 
the  equations  of  motion  suggests  that  we  put 


(5.18) 
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11  _ 11 

..1  n..£ 


nl  = 

n2 

= 0 

9 

2 

2 

= 0 

1 

1 

r 

C 

= n.l 

- n.3 

9 

P-2  = 

P-3  = 

P, 

12 

21 

= 0 

12 

21 

0 

rr 

= IT 

/V 

9 

£ = 

£ = 

22 

22 

11 

11 

tt  i = 
• *1 

"..2 

= 0 

9 

p..l 

II 

V 

k 

= 

i = p:3  = ° > 


22  22 

P n = P o = 0 

• *_L  • • d 


(5.19) 


Using  the  assumptions  (5.19) > the  functions  (5.5)  can  be  written  in  terms  of 
their  nonvanishing  components  as 

11  2 2 12 
2,  n3S3  » £ ".l^L  * £ ~ n.2^2  ’ n.l  n-2  ’ 


2 2 


P " P-i£l  5 £ " p.2%  5 r-l  ^ i'«2 


1 2 
P , = P 


(5.20) 


and 


11 

22 

11 

11 

22 

n = 

n 

= n -3% 

, TT  >3 

n .3  , 

11 

22 

11 

11 

22 

P 

r*J 

£ 

= P 

•3~3 

• p *3 

" P -3  • 

iform 

"pressure"  gradient  in 

the  fluid, 

(5-21) 


as  well  as  dp^/d^  are  independent  of  5,  i.e., 

3n  dp 

—1=0  — = 0 

35  ’ 35 

Also,  an  examination  of  the  equations  of  motion  (2.24)  and  (2.25),  together 
with  (5.19)  to  (5.21),  easily  reveals  that 


(5.22) 


so  that 


fl  - f2  - 0 * 


<12  - ‘13  ■ ^1  ■ & = 0 - 


11  11  22  22 

t . = i _ = r . = A 0 = 0 

••1  .*2  «.l  .*2 


(5-23) 


— r 
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In  view  of  the  assumed  forms  of  n,n°  and  na^  in  (5-19)  and  (5-20),  it  is 
at  once  apparent  that  the  conservation  equation  (2.26)  is  identically  satisfied. 
Since  the  nonvanishing  components  of  the  response  functions  (5«7)  have  the  forms 
(5-17)  and  are  therefore  independent  of  5,  it  follows  from  (2.17),  (2.20)  and 
(5.18)  that  at  time  t 

an 

+ Xf  3 = 0 . (5.27) 

Recalling  (5.20)^  and  (5.25)-^,  from  (5.22)^  we  obtain  aq/a5  = 0,  which  upon 
integration  yields 
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q = A1(t)?  + A2(t) 


(5.28) 


and  hence 

n3  =-a2(A1(t)^  + A2(t))  , \f3  = a^t)  , (5-29) 

where  the  coefficients  and  A 2 are  functions  of  t only.  Consider  next  the 

conservation  equations  (2.24)  and  recall  that  qa=  j£°  specified  by  (5.18)-. 

~ 2 

Since  tt0  and  p**  are  independent  of  F and  since  tt®  may  be  a function  of  F,  the 

equation  for  a-  1 in  (2.24)  after  using  (5.20)  and  (5*24)1  can  be  written  as 

v 1 

-1  -i  At 

-aT+w.i-".i*".i  • <5'30) 

Now  with  the  use  of  (5*25)g  ^ and  (5»22)2  and  by  argument  similar  to  that  which 
resulted  in  (5*28)  and  (5*29)  we  obtain 


= - 2av(t-To)ql;L  = B1(t)5  + B2(t)  , 


(5.31) 


U^1  = V*)  “a[Al(t)?"A2(t)] 


(5.32) 


Since  the  left-hand  side  of  (5. 31)  vanishes  at  t = tq  for  all  F,,  it  follows 
that  B1(tq)  = Bg ( To ) - 0. 

Again,  in  a similar  fashion,  from  (2.25),  ( 5 -23 ) > (5-24),  (5*26)  and 
(5.31)  we  obtain 

T (t-To>p“3  ‘ Bi(t)5*B2<t)  • 

2v  {+  t ,MU  2v  (+  t W11  R I11  - .22 

T (t_  o;U**3  T (t  0)  .«3‘ B].(t)  » *..3  ” *..3  • 


(5.33) 


The  expressions  (5.29)^  2,  (5.3J-)1»  ( 5 • 32 ) and  (5.34)  involve  the  four  undeter- 
mined functions  A^A^B^B,,.  We  determine  these  in  the  next  sections  by  an 


appeal  to  corresponding  results  in  the  three-dimensional  theory. 
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Before  closing  this  section,  we  observe  that  consistent  with  (5.22).. 

1>2 

and  corresponding  to  uniform  "pressure"  gradient  in  the  fluid,  the  difference 


in  the  values  of  n„  at  any  two  different  sections  of  the  pipe  is 
o 


An  = n - n = - a A A?  , 

A5=  ?2-5i 


(5.35) 


while  the  difference  AP#1  is  given  by 


4 * p!i[  -^1  , “ V5  • ■ VTo)45  ■ 0 

1 5=5?  5-5, 


(5.36) 


6.  Determination  of  the  unknown  coefficients  and  the  constitutive  response 
functions  in  the  solution  of  section 

This  section  is  concerned  with  the  determination  of  the  unknown  coefficients 

which  occur  in  (5.29)-^  2 anc*  (5*31)  to  ( 5 • 3^+ ) » as  well  as  the  relation- 
ships between  the  constitutive  response  functions  (5*7)  and  the  corresponding 
results  calculated  from  the  three-dimensional  solution  of  Poiseuille  flow.  We 
recall  that  the  final  results  in  section  5 make  use  of  the  fact  that  the  gradient 
of  the  axial  force  n^  along  the  center  line  of  the  pipe  is  uniform  and  this  cor- 
responds to  a uniform  "pressure"  gradient  in  the  fluid.  A result  of  this  kind, 
if  necessary,  enables  one  to  identify  the  coefficients  A^  and  B^  in  (5*35)  and 
(5.36)  experimentally  after  an  appeal  to  the  expression  for  the  resultants  n and 
p°  in  (A20).  Here,  however,  we  discuss  the  determination  of  the  unknown  functions 
A^,A2,B^,B2  through  a comparison  of  certain  expressions  in  section  5 with  cor- 
responding results  in  the  three-dimensional  solution  of  non-Newtonian  Poiseuille 
* 

flow. 

Preliminary  to  our  main  objective,  we  need  to  recall  certain  expressions 

and  results  from  the  three-dimensional  theory.  Let  (e  ,e  ,e_)  denote  the  unit 

''O 

base  vectors  of  cylindrical  polar  coordinates  (r,9,z)  and  recall  the  relationships 


e = cos  0 e,  + sin  0 e.  , e„  = - sin  0 en  + cos  0 e„  , 

'XL  'XL  '*'0  'Xl.  'XL 


(6.1) 


e,  = cos  0 e - sin  0 e.  , e_  = sin  0 e + cos  0 e„  , 

~1  ~r  ~0  ~9 

where  as  in  section  5 the  unit  base  vectors  (e..  ,e0)  ar^  us<d  in  place  of  (i,j). 

'XL  rJd  /x  ^ 

For  the  axisymmetric  problem  under  discussion,  the  Cauchy  stress  tensor  referred 
to  (e  can  be  written  as  [see  (A9)  of  Appendix  A] 

^ ^3  , 


An  account  of  non-Newtonian  Poiseuille  flow  may  be  found  in  [16,  Sec.  19]. 
The  solution  to  such  viscometric  flows  within  the  scope  of  the  theory  of 
viscoelastic  fluids  was  originally  given  by  Rivlin  [17]. 
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1 


T = T e 8e  + T (e  ®e  +e  ®e  ) 
~ <rr>~r  ~r  <r3>  ~r  ~3  ~3  ~r 


+ T e ®e  +T  e„®e_  , 

<99>~Q  ~9  <33>~3  ~3 


(6.2) 


where  ^<rr>’T<rz>>,1'<QQ>>T<zz>  are  nonvaniEhinS  physical  components  of  the 
stress  tensor  T in  cylindrical  polar  coordinates.  For  the  three-dimensional 
Poiseuille  flow,  these  stresses  are  given  by  [l6,  p.  20]: 

T<r3>='T  (K  } =-S  * S = 2 f ’ 


^r>-T<99>=^l(K)  ’ 


(6.3) 


T<33>  ‘ T<99>  a2 


* * 

= a0(K  ) , 


-X-  #■  % 

where  t is  the  shear  stress  function,  and  are  the  normal  stress  functions, 
Ak  * 

the  functions  am  are  related  to  am,  (m=l,2),  by 


*.  *,  *r  *-l.  A*.  . 

am(K  ) = am[T  (S  )]  = am(S  ) , (m-1,2)  , 

the  rate  of  shear  K is  defined  by  (4.31),  the  scalar  q can  be  computed  once 

* 

the  stresses  are  known  on  a cross-section  z = constant  and  f , called  the 
specific  driving  force,  is  a constant  given  by 


* 

F = 2tt 


na  (zg-z.^ 

0TH*2'T<33>k’rdr 


* * * 


(6.4) 


The  three-dimensional  viscometric  functions  t and  the  scalars  q ,f 

correspond,  respectively,  to  TjO-^jOg  and  q,f  of  [16,  Sec.  10].  For  further 

* 

properties  of  the3e  functions,  including  the  existence  of  the  inverse  of  t , 
see  [16,  Sec.  11].  In  particular,  we  note  that  by  virtue  of  a normalization  of 


the  constitutive  response  functionals , the  viscometric  functions  all  vanish 


X \ 

for  v = 0 (and  hence  for  K =0). 

Using  (A5)n  along  with  (Al8)  and  (4.28),  the  base  vectors  g and  their 

l,o 

reciprocals  g1  are  calculated  to  be 

1 

g.  = ae.  )9  e.  = a(cos  0e  - sin  0 e ) - 2v(t-t  ) -cos0e  , 

SI  o ~o  ~0  o a 


& = a^2 


2v(t-t  )02e_  = a(sin  9 e + cos  9 e ) - 2v(t-t  ) 7 s^n  ® eo  > (6.5) 

o a 


I 7 V 

o 


§>3  = ~3 


g1  = - e.  = - (cos  0 e - sin  9 e ) , g2  = 7 e = 7 (sin  0 e + cos  0 e ) , 

2.  a ~1  a ~0  2-  a a ~r  ~9 

S,3  = T (T-To)[0lf.l+  + % = T (T'To)r2r  + % • 

a 

Then,  by  (A10)  , the  vectors  T1  at  time  t may  be  expressed  as 

t t - 008  9 T<rr>Sr  ‘ 8in  9 T<99>5*  + 008  8 T<r3>t3  ’ 

k i ■ Bln  9 T<rr>Sr  + 808  9 T<99>d>  + sin  9 T<r3A  - 
T3  - t2v(t-T0)rT<[.r>+B2T<r3>]er  + [2v(t-T0)rT<r3>+a2T<33>]e3  . 

Clearly  with  the  use  of  (6.1).  _,  the  above  results  may  also  be  expressed 

1 yd 

relative  to  the  basis  (e^e^e^). 

Now  with  the  use  of  (6.7),  the  definitions  (A20)  and  following  a lengthy 
calculation,  the  resultant  forces  at  time  t are: 

. 2,  * *\  f ,eL2+r2  A*  _ A*\  , 

n = n{a  (f  5 + q ) - J (— jT-  0l-2ra2)dr}e3  , 


(6.6) 


(6.7) 


.1  1 


2 2 


H = n.l^L  ’ H _ tt.2^2  * n-l  n-2  9 

1 r P8  1 A*  * . . f*  *. 

tt.i  = na[-J  -Ojdr  + ^f  +v(t-TQ)— +q]  , 


(6.8) 


(6.9) 


1 1 


2 2 


P = P ,6,  » P = P.0ec 


2 a 4 4 . * 

1 2 Tta  ,*  a ^ T r -a  A*  A ^ v f 

3.rp.2  = ?H  Jn —r- aidr  + ^f +v(t_To)  t q ]]  » 

8 vJ 


(6.10) 


11  22  tt  2 * 12  21  . 

n =rr  =-r  a f e,  , tt  =tt  =0  , 

4 ivj  rv  ^ 


(6.11) 


su - e22 - i tfi-^ M)-r^y £%  , 

a 0 

12  21  _ 
p = p = 0 . 


(6.12) 


Also,  from  (A24),  (6.7) j the  expressions  for  Xf,XXa,XX° ^ at  time  t are: 


2 * 

\f  = - TTa  f e_  , 


(6.13) 


11  p p IP  r 1 A*  -X-  -X- 

XX  = XX. ^ , XX  = XX. 2^2  , XX. -L  = XX. 2 = TTa[-J  - c^dr  + gf  +q  ] , (6.l4) 


XX11  = XX22  = - 2§-  f*e  , XX12  = XX21  = o 

sss  (V  ^ AOJ  /V  n>J 


(6.15) 


It  is  worth  observing  that  while  the  motion  is  steady  both  in  the  solution 

via  the  three-dimensional  theory  and  the  direct  theory,  and  while  the  three- 

•* 

dimensional  stress  field  T(x,y,z,t)  is  also  steady,  the  expressions  for  1“  in 
(6.7)  and  for  the  resultants  n,p°,  etc.,  obtained  from  (A20)  involve  the  time  t 
explicitly.  This  is  due  to  the  fact  that  the  vector  in  (6.7)  and  hence  the 
resultants  n,pa,  etc.,  are  calculated  relative  to  the  convected  coordinate 
surface  g = constant  which  varies  with  z and  t,  as  can  readily  be  seen  from 
Fig.  2 or  Eq.  (4.30).  In  this  connection,  it  may  be  recalled  that  the  base 
vector  g^  which  is  normal  to  the  surface  f = constant  varies  with  time  but  is 
independent  of  both  g and  z. 

Now  the  coefficients  A^  and  B.^  in  the  solution  of  section  5 may  be 


35 


determined  by  requiring  that  the  difference  resultants  An^  and  kp  ^ in 

(5-35)  and  (5.36)  be  the  same  as  the  corresponding  expressions  calculated 

from  (6.8^  and  (6.10)  ,.  In  this  way,  we  obtain 

3 


* na  / \ * 

Ax  = - rtf  , Bx  = — v(t-TQ)f 


(6.16) 


By  (5.29),  (6.16)^  and  the  fact  that  n (0)=0,  in  the  solution  of  section  5>  we 
have  n^(0) = n^(0) = -a  k^.  Comparison  of  this  with  the  corresponding  result 
calculated  from  the  three-dimensional  solution  yiexds 


A2  = - rrq 


(6.17) 


It  remains  to  identify  the  coefficient  B^.  This  can  be  effected  with  the  use 
of  (5.31)  and  (5-33)  and  by  comparing  the  results  for  p^  and 


B2=fv(t-To)q  . 


(6.18) 


An  examination  of  the  various  results  also  reveals  that  the  Lagrange  multiplier 

n pp 

q =q  which  occurs  in  (5.22)  is  not  an  independent  quantity  and  is  related 
to  the  multiplier  q in  (5.25)  by 


11  1 
q = i q 


(6.19) 


Having  made  the  identifications  (6.16)  to  (6.18),  the  constitutive 
response  functions  (5.7)  may  be  also  identified  from  the  corresponding 
expressions  in  the  three-dimensional  solution.  We  omit  details  here  but  note, 


in  particular,  the  expression  for  n^  given  by 


A f rS^+r2  A*  „ 

n3  =_n  0l"2ro2]dr  ’ 


(6.20) 


which  clearly  relates  the  effect  of  "normal  force"  to  the  "normal  stress" 


effect. 
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Appendix  A 


We  record  in  this  appendix  certain  details  of  an  approximation  procedure 
[6,7,8]  whereby  the  basic  equations  governing  the  motion  of  jets  (or  rods)  can 
be  derived  by  integration  from  the  three-dimensional  equations  of  classical 
continuum  mechanics.  Although  the  integrated  equations  have  been  obtained 
previously  [6,7,3]  for  any  number  of  directors,  we  record  here  only  those  which 
correspond  to  a directed  fluid  jet  with  five  directors.  We  also  record  the  con- 
straint conditions  appropriate  for  the  incompressible  Poiseuille  flow  discussed 
in  the  present  paper. 


First  we  define  what  is  meant  by  a jet-like  body.  For  this  purpose, 

consider  a finite  three-dimensional  body  © in  a Euclidean  3-space,  and  let 

convected  coordinates  Q1  (i= 1,2,3)  be  assigned  to  each  particle  (or  material 

+ * 

point)  of  ffl.  Further,  let  r be  the  position  vector,  from  a fixed  origin,  of 

a typical  particle  of  © in  the  present  configuration  at  time  t.  Then,  a motion 

A* 

of  the  (three-dimensional)  body  is  defined  by  a vector-valued  function  r which 

* S 

assigns  position  r to  each  particle  of  ffl  at  each  instant  of  time,  i.e., 


* A*,„l  „2  ft3  . \ 

r = r (0  ,0  ,0  ,t) 


(Al) 


We  assume  that  the  vector  function  r --a  1-parameter  family  of  configurations 
with  t as  the  real  parameter  — is  sufficiently  smooth  in  the  sense  that  it  is 
differentiable  with  respect  to  and  t as  many  times  as  required.  It  is 

•3 

convenient  to  set  0 - | and  adopt  the  notation 


91  = (9°, 5)  > 93  = 5 


(A2) 


We  shall  be  concerned  here  with  material  curves  (not  necessarily  straight  .lines) 


+The  use  of  an  asterisk  attached  to  various  symbols  is  for  later  convenience. 
The  corresponding  symbols  without  the  asterisks  are  reserved  for  different 
definitions  or  designations  to  be  introduced  later. 

^Recall  that  when  the  particles  of  a continuum  are  referred  to  a convected 
coordinate  system,  the  numerical  values  of  the  coordinates  associated  with 
each  particle  remain  the  same  for  all  time. 
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in  ® defined  by  the  equations  0 = 9 (^),  (a =1,2);  the  equation  resulting  from 
(l.l)  with  0O=8°(^)  represents  the  parametric  form  of  this  material  curve  in 
the  current  configuration  and  describes  a 1-parameter  family  of  curves  in 
space,  each  of  which  we  assume  to  be  smooth  and  nonintersecting.  We  designate 
the  space  curve  9®=0  in  the  current  configuration  by  c.  A point  of  thi6  curve 
is  specified  by  the  position  vector  r,  relative  to  the  same  fixed  origin  to 

r*J 

* 

which  r is  referred,  where 


A . A*,  . 

r = r(|,t)  = r (0,0, §,t)  , 


(A3) 


with  5 belonging  to  a finite  interval  [ ^ , ^ ] . Let  a^a^  anci  a3  den0‘te  the 
unit  principal  normal,  the  unit  binormal  and  the  tangent  vector,  respectively, 
to  the  curve  c.  At  each  point  of  c imagine  material  filaments  lying  in  the 
normal  plane,  i.e.,  the  plane  perpendicular  to  a_,  and  forming  the  normal  cross- 
section  0.^.  The  surface  swept  out  by  the  closed  boundary  curve  dfln  of  (2 n is 
called  the  lateral  surface.  Such  a three-dimensional  body  is  called  jet-like 
if  the  dimensions  in  the  plane  of  the  normal  cross-section  are  small  compared 
to  some  characteristic  dimension  L(c)  of  c (see  Fig.  l),  e.g.,  its  local  radius 
of  curvature  l/k,  or  the  length  of  c in  the  case  of  a straight  curve.  A jet- 
like body  is  said  to  be  slender  if  the  largest  dimension  of  dn  is  much  smaller 
than  L(c).  If  dn  is  independent  of  the  body  is  said  to  be  of  uniform  cross- 
section,  otherwise  of  variable  cross-section.  Let  the  (three-dimensional)  jet- 
like body  in  sene  neighborhood  of  c be  boundary  by  material  surfaces  5 = 5^, 

5 =§2  (indicated  in  Fig.  1)  and  a material  surface  of  the  form 

F(e1,e2,?)  = 0 , (A4) 

which  is  chosen  such  that  5 = constant  are  curved  sections  of  the  body  bounded 


The  normal  cross-section  of  a jet  is  a portion  of  the  normal  plane  to  the 
curve  c,  i.e.,  the  intersection  of  the  body  and  the  normal  plane. 


4o 


r 


i 


by  closed  curves  on  this  surface  with  c lying  on  or  within  (A4).  In  the 
development  of  a general  theory,  it  is  preferable  to  leave  unspecified  the 
choice  of  the  relation  of  the  curve  c to  the  boundary  surface  (a4).  In  special 
cases  or  in  specific  applications,  however,  it  is  necessary  to  fix  the  relation 
of  c to  the  surface  (A4). 


(A5) 

(A6) 


* 

and  further  assume  that 

e5  = > 0 * (A7) 

In  (A5)-(a6),  and  g1  are  the  covariant  and  the  contravariant  base  vectors  at 

time  t,  respectively,  g. . is  the  metric  tensor,  g^  is  its  conjugate,  6^  is  the 

d 

Kronecker  symbol  in  3-space  and  dv  the  volume  element  in  the  present  configuration. 

* 

The  velocity  vector  v of  a particle  of  the  three-dimensional  body  in  the 
present  configuration  is  defined  by 

* •*  , 

v = r , (A8) 


where  a superposed  dot  denotes  material  time  differentiation  with  respect  to  t 

holding  Q'1'  fixed.  The  stress  vector  t across  a surface  in  the  present  configura 

# 

tion  with  outward  unit  normal  v is  given  by 


++For  most  purposes, we  could  assume  a less  general  form  for  the  lateral  bounding 
surface  of  the  body  and  write  F(0^,9^)  = 0 instead  of  (A4). 

The  choice  of  positive  sign  in  (A7)  is  for  definiteness.  x Alternatively,  for 
physically  possible  motions  we  only  need  to  assume  that  g2  ^ 0 with  the  under- 
standing that  in  any  given  motion  is  either  >0  or<0.  The  condition 
(A7)  also  requires  that  0*  be  a right -handed  coordinate  system. 


where 


* i\r  ---1  ik 

t = Vig2T  = viT  & , T = g.®g2T  =T  g.®^  , 


t1  = gMV  = gK]gJ 

~ "0  J~ 


* i 
= v.g 


g1  * TgJ 


(A10) 


T is  the  symmetric  Cauchy  stress  tensor  and  ® denotes  the  tensor  product  of  two 
vectors.  In  terms  of  quantities  defined  in  (A5)-(A10),  the  local  field  equations 
which  follow  from  the  integral  forms  of  the  three-dimensional  conservation  laws 
for  mass,  linear  momentum  and  moment  of  momentum,  respectively,  are 


* i 

P g2  = 0 , 


(All) 


T + p f g2  = p g2v  , g.  XT  =0  , 

~ jl  /^-tL  'X'  r*t  * 


where  p is  the  three-dimensional  mass  density,  f is  the  body  force  field  per 

unit  mass  and  a comma  denotes  partial  differentiation  with  respect  to  e1. 

Let  a (not  necessarily  the  same  as  the  normal  cross-section  Ctn  defined 

above)  denote  the  curved  section  of  the  surface  g = constant  bounded  by  da,  i.e., 

a closed  curve  which  is  the  intersection  of  the  surface  g = constant  and  the 

surface  (,AU).  ttirther,  let  the  points  g^  and  g^,  with  g^<g2,  ^orm  endpoints 

of  a segment  of  the  curve  (A3)  which  we  denote  by  P,  and  designate  by  ( 1^  and 

the  particular  sections  associated  with  g^  and  g^,  respectively.  Now  consider 

* 

an  arbitrary  part  of  p of  the  three-dimensional  region  occupied  by  the  body 

* * X- 

such  that:  (i)  P contains  p;  and  (ii)  the  boundary  dP  of  P consists  of  the 

sections  0^  and  0-,  and  a portion  of  the  surface  (A4)  bounded  at  each  end  by 

* 

dax  and  da2-  A body  so  described  is  called  a jet-like  body  and  the  part  P 
forms  a portion  of  such  a body. 


The  mass  m of  a portion  of  the  jet-like  body  is  given  by 


t.  *,  p * 

(**  ) = Jp*  P dv 


* r * 1 i p 

dv  = * p g5de  de  dg  . 

Jp 


(A12) 


where  p is  the  mass  density  of  the  (three-dimensional)  continuum  and  dv  is  the 
element  of  volume  in  the  present  configuration  at  time  t.  In  terms  of  the 
segment  P of  the  material  curve  (A3),  the  same  mass  has  the  alternative 
representation 


m 


f [f  p g2d01d02]d|  = I"  p(a  )2d$  = [ pds  = m(p)  , 

Jp  Jp  Jp 


where  the  density  p per  unit  length  of  the  curve  (A3)  is  defined  by 


(A13) 


x = p(a33)2  = J p g2deW 

In  view  of  (All)^,  we  note  that 


(A14) 


'Q 


X = 0 


The  curve  (A3)  is  fixed  in  the  jet-like  body  by  the  condition  [6] 


(A15) 


f p*g2e°de1 
Ja 


d0  = o . 


(Al6) 


(A17) 


As  in  the  paper  of  Green  and  Naghdi  [7],  we  assume  that  for  the  jet-like 

* . . 

body  described  in  this  appendix,  the  position  vector  r in  (Al)  can  be  repre- 
sented  by  the  expansion 

*12  x ®i  °2  aN 

* (0  ,0  ,?,t)  = r + e 0 ■L0  *...0  “d  „ , 

~ ~ N ~ai<V'*aN 

where  r and  d are  vector  functions  of  f,t.  In  (A17),  the  vectors 

d are  completely  symmetric  in  the  indices  o, jou, • . . ,aM,  the  summation 

• • • OJjy  -L  c-  N 

is  over  all  values  of  <*^>0, » • • • ,Ojj  and  N = l,2,3,...  We  assume  that  (A17)  may 
be  differentiated  as  many  times  as  required  with  respect  to  any  of  its  variables. 

In  the  rest  of  this  appendix  we  consider  a special  case  of  (A17)  and  restrict 
attention  to  the  approximation 


‘(01t02»5.t)  = r(|,t)  + 0^(5, t)  + 0®0Pd^(s,t)  • 


(A18) 


T*'  - -s*..*.* • 
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Using  this  assumption  in  (Al),  (A5)  and  (A8),  we  obtain 


* X a 0 ‘ ’ 

v = v + ©*w  + 0 0Kw  „ , v = r , w = d , 

- - ~o  ~a&  ~ - -a  -a 

g,  = d +2Qed_  > go  = ao  + eV+9aeV.0  > 


-a  -a 


-aB  ’ -3  ~3  -o 


~o0 


w _ = d „ , 

-oB  — o0 


a = r 


(A19) 


where  a prime  denotes  partial  differentiation  with  respect  to  5.  We  also 
recall  the  following  definitions  for  the  resultants  nj-n^p®,!/*0,?®0: 


/'■»>  /v 


n = f T3d01d02  , tt®  = f TQfd01d02  , p®  = [ ^©“d©1. 

~ Ja  ~ JQ  ~ Ja  ~ 


d02  , 


tt®0  = a,®0**0"  , .*«[  Ta0pd01d02  , p^=f  T30a9Pd01d02  . 

~ Ja  ~ ~ Ja  ~ 

The  equations  of  motion  in  terms  of  the  resultants  (A20)  are  obtained  by 

suitable  integration  of  (All),,  over  a section  d and  are  given  by  (for  details 

>3 

see  [6,7]): 

Bn  dpa  o 

, -£-♦  Xi*=  n^Xy*9^  , 

~ ~ <35  ~ ~ 

3D  0(0 

— . . .ckB  aB  , , a0v6' 

“ST + ~ ~ + ’ 


(A20) 


r,  X n + d x tt®  + d X tt®0  + d/  X P®  + d/  x P®0  = 0 , 

provided  that 

\f  = f p*g^f  "d01d02  + f [(d©2!1  - doV)  + T3(\2d01  - X^'d©2)  3 , 

\jia  = f 0ap  Vf  Wd©2  + f ©*[  (d02Ti  - deV ) + T3  ( X2d0X  - xW ) ] , 

~ Ja  Jaa  - - - 

Xi®0  = f ©“©^pV^f  Wd©2  + r e“©P[ (d©2!1  - d©1!2)  + T3(x2d01  - X^©2) ] 

Ja  ~ Jaa  - - - 

and 


(A21) 

(A22) 

(A23) 


' a 


(A24) 


xy°*  - f (*sW«W  , x,*'1* 
Ja 
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f p*g^0“©P0Y06d01d02  , 

Ja 


(A25) 


T 


where  \“  = \ • ga  and  X is  a vector  tangential  to  the  boundary  surface  (A4). 

/x  rs* 

If  we  adopt  the  approximation  (Al8)  and  identify  the  vectors  d ,d  and 

'■-a  ~ofp 

the  position  vector  r in  (Al8)  with  the  directors  in  ( 2.1 and  the  position 
vector  (2.1)^  of  the  curve  c,  then  the  development  of  this  appendix  and  the 
corresponding  results  given  in  section  2 are  formally  equivalent-  In  particular , 
comparison  of  the  equations  (A15)  and  (A2l)  to  (A23)  with  those  in  (2.32) 
reveals  a 1-1  correspondence  between  the  two  systems  of  equations  provided  we 
identify  the  expressions  (A24)  to  (A25),  respectively,  with  the  assigned  fields 
and  the  inertia  coefficients  in  the  theory  of  a directed  curve  with  five 
directors  discussed  at  the  end  of  section  2. 

Before  closing  this  appendix,  we  discuss  an  appropriate  constraint  con- 
dition arising  from  incompressibility  when  the  position  vector  is  approximated 
by  (A18).  For  a three-dimensional  incompressible  medium,  the  mass  density  p 
is  constant  and  by  (All)^  the  condition  of  incompressibility  at  time  T^t  is 

£ l&Lfcfel  - 0 > <«= 6> 

1 

where  g2  is  defined  by  (A7)^.  If  we  confine  attention  to  the  special  case  in 
which  the  position  vector  is  specified  by  the  approximation  (Al8),  then  with  the 
use  of  (A5)1,  we  find 

e5  = [iifes1  * t2eXV[i)£Yvfe]1 

+ 9>'9Vt[iifcixv]+2«,‘v[did)iYd;] 

+ ‘•tiuSvaJa1!  t29°'9'l9Vt2[2ai3i2i:1 
" eP“[iP4«,ii:v1>  + 1*9“9S9’'9Y[5«i2e22l:vI  ■ (A27) 


where  e^v  and  the  operator  ( )'  are  respectively  defined  by  (3*5)  and  (5«10). 
After  substituting  the  above  result  in  (A26)  and  following  a routine  calculation 
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45. 


we  obtain  the  following  conditions 


dT  C£i^3J  = 0 ’ (A28> 

h +2e^v  £ fixivvS.3]  * 0 ■ <*»> 

^[4i&V+e"v^l[^~Tiv^I  + [Vvvi;it2tiiiAvi3n‘0  • (A30) 

+ + [ + tVfSa^1  + [2apVe’ 

+ tV*&n'°  • (A31) 

ePCT  -f-  {[d  d^  d'  ]+td  da  d'  ] + [dQ  d,,  d'  ] 
dT  ~'Yp~T]a~Qrf3  ~Vp~^o~o?T|  ^p~Hor<r/ 

+ [d  d d']  + [d  d d'  1 
~ap~'Ya~&7]  ~op^a~YT] 


+ [d  d d'  ]}  = 0 . 

~ap~T|a~$Y  J 


(A32) 


On  the  basis  of  the  1-1  correspondence  noted  in  the  preceding  paragraph,  we 
may  employ  the  five  conditions  (A28)-(A32)  as  constraint  conditions  for  the 
special  motion  of  tne  theory  of  a Cosserat  curve  of  an  incompressible  fluid 
discussed  in  section  5.' 


Figure  1:  A jet-like  body  in  the  present  configuration  showing  the  curve  c 
with  position  vector  r (defined  in  Appendix  A)  as  the  line  of 
centroids  and  the  end  normal  planes  s0  6hown  are 

the  unit  principal  normal  a^,  the  unit  binormal  and  the  tangent 
vector  a«  to  the  curve  c. 


d2(C|,T0)=aj 


d2^,r0)=aj< 


d2(e2,r0)  = aj( 


£z=al 


r.  (f  o , r) 


ai|(T)  = d22(r) 


Figure  2.  A special  Poiseuille  flow  in  a circular  pipe:  Fig.  2(a)  depicts  a 

reference  configuration  at  time  tq;  Fig.  2(b)  depicts  the  corresponding 
flow  configuration  at  time  t;  and  Fig.  2(c),  along  with  Fig.  2(b), 
indicate  the  locations  A'  and  B'  of  particles  on  the  center  line  of 
the  pipe  which  at  time  tc  were  at  A and  B,  respectively. 


r 
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Figure  3:  A sketch  of  an  enlarged  portion  of  Fig.  2(c)  for  the  upper-half  of 

a cross-section  of  the  flow  in  the  x-z  plane  showing  the  position 
# 

vector  r (t)  of  a fluid  particle  at  time  t. 
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